Neutron scattering study of novel magnetic order in Nao.5Co02 
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We report polarized and unpolarized neutron scattering measurements of the magnetic order in 
single crystals of Nao.5Co02. Our data indicate that below Tjv = 88 K the spins form a novel 
antiferromagnetic pattern within the C0O2 planes, consisting of alternating rows of ordered and 
non-ordered Co ions. The domains of magnetic order are closely coupled to the domains of Na 
ion order, consistent with such a two-fold symmetric spin arrangement. Magnetoresistance and 
anisotropic susceptibility measurements further support this model for the electronic ground state. 
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PACS numbers: 71.27.+a, 74.70.-b, 75.25.-(-z 

The layered cobaltates Na2;Co02 have attracted much 
attention due to their unusual thermoelectric proper- 
ties T,"^ and because of the recent discovery of supercon- 
ductivity in the hydrated composition ^ . The Co ions in 
these compounds form a hexagonal lattice, and the aver- 
age valence can be changed by varying the Na concentra- 
tion X. The Co^+ ions can formally be regarded as mag- 
netic ions with S* = on a frustrated low-dimensional 
lattice. It is believed that strong electronic correlations 
play an important role in the physics, and, as a function 
of the Na concentration, a rich phase diagram has been 
reported For X < ^, the material is a paramagnetic 
metal, while for a; > ^ an unusual "Curie- Weiss metal- 
lic" phase is observed. At a; = i, a unique state is real- 
ized: the Na ions are chemically ordered to form zig-zag 
chains in an orthorhombic superstructure with two-fold 
symmetry |^ 0], and it has been speculated that this 
leads to Co^'^+VCo'^'^~* charge order within each C0O2 
layer H HH H The ground state is believed to 

be a magnetically ordered insulator; resistivity, Hall co- 
efficient, thermal transport 0|, and angular magnetore- 
sistance oscillation measurements are all consistent 
with such a two-fold symmetric ground state. 

In order to further investigate the electronic ground 
state of the half-doped C0O2 plane, we have performed 
neutron scattering, susceptibility, and transport studies 
on single crystal samples. The combination of polarized 
and unpolarized neutron scattering data allow us to de- 
termine the ordered spin direction and arrangement. We 
find that the ground state is well described by an or- 
dered array of stripes of antiferromagnetic spins inter- 
leaved with stripes of non-ordered Co ions. Our trans- 
port and anisotropic susceptibility data further support 
such a model. 

The samples used in this study were prepared by 
electrochemically deintercalating a floating-zone grown 
Naa;Co02 crystal to yield a final concentration oi x — 
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FIG. 1: Bulk physical properties measured using Nao.5Co02 
single crystals, (a) In-plane resistivity pab- (b) Fractional 
change of pab measured in a magnetic field of 14 T with -ff || c 
and H \\ ah. (c) Anisotropic magnetic susceptibility in a field 
of 1 T. The vertical dashed lines denote the temperatures of 
T = 51 K and Tjv = 88 K. 

0.5 [1^ Il3| . Figure 1 shows measurements of the 
bulk properties of the resulting Nao.5Co02 crystals. In 
Fig. 1(a), we plot the in-plane resistivity {pab) measured 
using the conventional four-probe method on a crystal 
with dimensions 1 x 0.5 x 0.05 mm'^ in zero applied field. 
The resistivity dramatically increases below 51 K, consis- 
tent with previously published data . The application 
of a 14 T magnetic field affects pab in different ways de- 
pending on the field orientation, as shown in Fig. 1(b). 
For H \\ ah, a negative magnetoresistance is observed be- 
low ~ 51 K, consistent with previous work by Balicas et 
al. 11]. The kink in the data near 50 K is caused by a 
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shift in the onset temperature of the resistivity upturn. 
A new observation is that, for H \\ c, the data exhibit 
a positive magnetoresistance below 88 K. As we discuss 
further below, this behavior coincides with the magnetic 
phase transition at T^v = 88 K and is consistent with 
our proposed ground state. The magnetic susceptibility 
(Fig. 1(c)) is anisotropic and exhibits kinks at T = 88 if 
and 51 K [HEl- 

Our neutron scattering experiments were conducted at 
the NIST Center for Neutron Research using the BT9 
and BT2 triple-axis spectrometers. The incident energy 
was fixed at 14.7 meV, and the coUimations were 40'-40'- 
sample-60'-open. Pyrolytic graphite filters were placed 
in the beam before the sample to reduce higher order 
neutron contamination. For the polarized measurements, 
Heusler crystals were used for the monochromator and 
analyzer, and polarization guide fields were created at 
the sample position using Helmholtz coils (a flipping ra- 
tio of ~ 21 was achieved). The crystal studied was cylin- 
drically shaped, with mass 2.41 g and a mosaic of 1 to 
2° (full- width half maximum), depending on the orienta- 
tion. The temperature was controlled using a ^He closed 
cycle refrigerator. We label reciprocal space peaks using 
the hexagonal PQ^/mmc space group, with low temper- 
ature lattice parameters a = 2.816 A, and c = 11.05 A. 
As mentioned previously, in Nao.5Co02, Na ion order- 
ing leads to an orthorhombic supercell (oo = V^a and 
bo — 2a) in real space, which produces nuclear superlat- 
tice reflections in reciprocal space 0. 

Figure 2 shows unpolarized neutron diffraction evi- 
dence for magnetic order below Tn — 88 K. At low tem- 
peratures, new peaks arise which are distinct from the Na 
superlattice reflections. In the (iJiJL)-zone, we observe 
such peaks at Q = (^ ^ odd). The temperature depen- 
dence of the peak intensities of the (| ^ 1) and (^ ^ 3) 
peaks are shown in Figure 2(a). Also plotted are the in- 
tegrated intensities of the (5 5 1) reflection which match 
the T-dependence of the peak intensity. Figures 2(b) 
and 2(c) show reciprocal-space scans through the (5 5 1) 
peak along the (HHO) and (OOL) directions, respectively, 
at T = 8 K and T = 120 K. The low temperature peak 
widths are resolution-limited, indicating long-range order 
(> 50 A). 

In order to confirm that these peaks originate from 
magnetic order, we performed polarized neutron mea- 
surements. This is the clearest way to distinguish mag- 
netic peaks from weak nuclear superlattice peaks. In gen- 
eral, the magnetic diffraction cross section is given by 
im) - E.,/ I (/I El e"3-^~ C/f"^^K)P, where t and / de- 
note the initial and final states of the system. Uj^'^^ ^ 
F{Q){S f\M±i ■ <T\Si) is the scattering amplitude for neu- 
tron spin state Si to Sf at atomic site I, where F{Q) is the 
magnetic form factor, AI± is the component of the mag- 
netic moment perpendicular to the scattering wavevector 
Q, and a is the spin operator. 



100 120 




-1.2 -1.0 



FIG. 2: (a) Temperature dependence of the intensities of the 
(1 1 1) (°) ^^'^ (5 1 (^) magnetic Bragg peaks mea- 
sured by neutron diffraction (a background was subtracted). 
The filled triangles (►) denote the integrated intensities ob- 
tained from 6-scans, and the filled circles (•) denote the in- 
tensities measured using polarized neutrons. The solid line 
corresponds to I oc (M>^ oc (Tjv - Tf^ with Tjv = 88 K and 
P = 0.28(1). Inset: the ratio of the intensity of (^ ^ 1) to that 
of (^ I 3). (b) (HHl) scan and (c) (| ^ L) scan through the 

magnetic Q = (| ^ 1) peak at T = 8 K (•) and T = 120 K 
(o). The solid lines are Gaussians fits, and the horizontal bars 
indicate the experimental resolution. The lower panels show 
polarized neutron diffraction with (d) P \\ Q and (e) P _L Q, 
for spin-fiip (•) and non-spin-flip (o) channels at T = 8 K. 



We have measured both the spin-flip and non-spin flip 
cross-sections with neutron polarizations both parallel 
and perpendicular to the wavevector Q. Figure 2(d) 
shows the scattering measured at T = 8 K with P \\ Q 
at the (i i 1) peak in the {HHL)-iomi. In this ge- 
ometry, all magnetic scattering will occur in the spin- 
flip channel, and nuclear scattering will be non-spin-flip. 
The data show that the scattering is entirely spin-flip. 
The small peak in the non-spin-flip channel is consis- 
tent with the background measured at high-temperature 
(not shown), taking into account the instrumental flip- 
ping ratio. The temperature dependence of this scat- 
tering (scaled to match the unpolarized data at low-T) 
is denoted by the filled circles in Figure 2(a). We see 
that the temperature dependence is identical. Hence, 
this proves that the scattering below 88 K at (^ ^ 1) is 
entirely magnetic. 

Neutron polarization analysis is also an extremely use- 
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FIG. 3: (a) Neutron diffraction peaks measured in the (HKO)- 
zone at T — 8 K. The three sets of data correspond to peaks 
from three different twin domains due to Na ion ordering: 
(i i 0) (•), (1 i 0) (o), and (i 1 0) (□) (the angle 9 was offset 
by 60° for two of the peaks). The weak peak at ~ 56° is due 
to a small crystallite, (b) Magnetic Bragg peaks measured in 
scattering zones tilted from the (_H'_R'0)-zone by 14.3° in order 
to reach the positions (i i 1) (•), (1 i 1) (o), and (i 1 1) 
(□). The solid lines correspond to Gaussian fits. 

ful method to deduce the direction of the ordered mo- 
ment. In Figure 2(e), we plot the scattering data mea- 
sured with P _L Q; that is, the {H H L)-zone of the sam- 
ple is horizontal and the neutron polarization direction 
is vertical. In this geometry, magnetic scattering can oc- 
cur in the non-spin-flip channel if the ordered moment 
has a component parallel to the polarization direction. 
A peak in the spin-flip channel would arise from ordered 
moments with a component in the (HHL) -plane. Here, 
the signal predominately occurs in the non-spin-flip chan- 
nel. Therefore, the direction of the ordered moments lies 
within the ab plane, and the ordered moment along the 
c-axis must be small (or zero) . The error bars on our data 
allow us to conclude that the size of the c-axis compo- 
nent which contributes to the scattering at (^ ^ 1) is less 
than ^ 0.04/is- To further specify the direction of the 
moments within the plane, the twin domain distribution 
must be taken into account, as we discuss below. The 
temperature dependence of the polarized data indicates 
that the c-axis component remains small at all measured 
temperatures below Tn = 88 K. 

The refinement of the magnetic structure is compli- 
cated by the presence of multiple twin domains (struc- 
tural and magnetic). These domains, in principle, can 
have different populations. To characterize the structural 
domains, diffraction measurements were taken with the 
sample aligned in the (7Ji^0)-zone. Since the Na ion or- 
der reduces the six-fold hexagonal symmetry to two-fold, 
there are three possible Na-order domains whose rela- 
tive orientations differ by a rotation of 60° about the 
c-axis. In Fig. 3(a), we show 0-scans through peaks at 

Q = h h 0)' ^^'^ (5 ^ Each peak is com- 

prised of scattering from the Na superlattice of two dif- 
ferent structural twin domains. The measured intensities 
of the peaks are proportional to the populations of the 
twin domains which contribute to that peak. (We note 



(a) Nuclear Bragg peaks 

rnucl rnucl 
^cxp ^calc 

(1 0) 561(12) 603 

(0 4) 624(27) 568 

(i i 0) 104(9) 96 

(I 2) 9.2(0.4) 10.0 



(b) Magnetic Bragg peaks 

rmagn rmcign 
-"cxp £calc_ 

(i ^ 1) 3.9(0.3) 3A~ 

(i I 3) 2.2(0.1) 2.3 

(i I 5) 0.6(0.3) 1.2 

(| 1) 8.1(4.4) 4.4 

(I 3) 3.9(1.6) 3.8 



TABLE I: Observed and calculated integrated intensities of 
selected nuclear and magnetic Bragg peaks at T = 8 K. The 
large error bars for the (^0 odd) magnetic peaks are due to 
the presence of nuclear superlattice peaks at the same posi- 
tions. For these peaks, the intensities at T = 120 K were 
subtracted from the intensities at T = 8 K. The calculated 
intensities are based on a model which includes the twin do- 
main distribution, an isotropic Co*^ magnetic form factor, 
and the instrumental resolution function. 

that the long axis of our cylindrically shaped sample lies 
in the scattering plane in this geometry. Hence, part 
of the intensity variation is related to neutron absorp- 
tion.) The Na order peaks are resolution limited in all 
three crystallographic directions, indicating that they are 
macroscopic in size. Also, the relative intensities of these 
superlattice peaks compared to the fundamental reflec- 
tions, shown in Table I, indicate that the Na ordering 
occurs in essentially 100% of our sample. 

In conjunction with characterizing the structural 
twins, we have also performed scans_ through the mag- 
netic positions: (| ^ 1), (1 ^ 1), and (^ 1 1). These scans, 
shown in Fig. 3(b), were performed by first aligning a 
(^ ^ 0) peak, then tilting the sample by 14.3° to bring 
the {5 5 1) position into the horizontal scattering plane. 
We find that the intensities of these magnetic peaks are 
proportional to the intensities of the corresponding Na 
order peaks. This suggests that each structural twin do- 
main corresponds to a magnetic twin domain. That is, 
the arrangement of ordered magnetic moments is closely 
coupled to the arrangement of the Na ions. This is a nat- 
ural expectation due to the symmetry of both arrange- 
ments and now has experimental support. 

The above information on the domain distribution 
and the moment direction can be combined with mea- 
surements of intensities of several nuclear and magnetic 
Bragg peaks to determine the magnetic structure. Our 
data are well described by the model shown in Figure 4. 
Within each C0O2 layer, for a single magnetic twin do- 
main, the Co ions are arranged in alternating rows (or 
stripes) of magnetically ordered and "non-ordered" ions. 
For the magnetic ions, the nearest neighbor spins are 
coupled antiferromagnetically, both along the row and 
between rows. Between C0O2 planes, the rows of mag- 
netic Co ions are stacked directly on top of each other, 
and the nearest neighbor interplane coupling is also an- 
tiferromagnetic. The directions of the ordered moments 
are parallel to the rows. The error bars on our data al- 
low us to put restrictions on the deviations of the ordered 



4 




FIG. 4: Model of the spin arrangement in Nao.5Co02. The 
two C0O2 layers in the magnetic unit cell are depicted, and, 
for clarity, only Co ions are shown. The solid circles repre- 
sent "non-ordered" Co ions, and the hollow circles represent 
magnetically ordered Co ions where the arrows indicate the 
directions of the magnetic moments. The shaded parallelo- 
gram outlines the magnetic unit cell. 

moments away from this direction to be less than 10° out 
of the ab plane. However, within the ab plane, our data 
are less restrictive, and the moments may be oriented up 
to 30° away from the stripe direction. The low temper- 
ature {T — 8 K) static magnetic moment is 0.25(2) 
per magnetic Co ion. The size of the observed moment is 
smaller than that expected for 5 = 1/2, which may sug- 
gest the presence of quantum fluctuations or deviations 
from a local moment picture. The order that we observe 
differs somewhat from the models proposed by Yokoi et 
al. [l^ , though we agree on the fundamental periodicity. 

The structure proposed in Figure 4 is the simplest 
model which is consistent with our data, assuming a local 
moment description of the magnetic order. It is consis- 
tent with having rows of Co^^ ions which are magnetic 
(in the S = 1/2 low spin state) alternating with rows 
of Co'^"'" ions (in the S = spin state) which are non- 
magnetic, as suggested previously However, we em- 
phasize that our neutron data do not directly probe the 
valence of the Co ions, only their magnetic moments. For 
example, the rows of non-ordered ions may correspond 
to Co ions with fluctuating moments which do not con- 
tribute to the long-range antiferromagnetism. The va- 
lence difference between the two distinct Co sites may, 
in fact, be small. Recent NQR studies suggest that 
unit-valence charge disproportionation does not occur 
at T/v 01 . It is possible that some degree of valence 
disproportionation exists already at temperatures above 
T/v = 88 K, and it is primarily the spin components which 
order at T/v 0|. A delocalized moment model may also 
describe the data, so long as the periodicity and spin 
orientation are the same as discussed. 

The anisotropic magnetic susceptibility measurements 



shown in Figure 1(c) are consistent with our model. A 
clear drop is seen at T/v = 88 K for _ff || ab, whereas 
no kink is visible for H \\ c, consistent with the obser- 
vation that there is essentially no c-axis component to 
the ordered moment below T/v. The susceptibility data 
also show an almost isotropic drop below T = 51 K, 
coincident with the upturn in Pab{T). A previous /iSR 
study suggested that a spin-reorientation transition oc- 
curs near this temperature We find that there is 
no indication of this transition in the magnetic neutron 
diffraction. Specifically, the ratio of intensities of the 
(^ ^ 1) and (^ ^ 3) Bragg peaks should change if there 
were a significant spin-reorientation below T = 51 K, 
but the inset of Fig. 2(a) shows that the ratio is con- 
stant. This is consistent with the polarized neutron re- 
sults and implies any spin reorientation must be small 
(< 10° out of the ab plane). Further, this suggests that 
the drop in susceptibility near T = 51 K is not caused 
by the spins which participate in the magnetic order be- 
low T/v = 88 K. Hence, this feature may be related to a 
gap developing in the excitation spectrum of the "non- 
ordered" Co ions. Finally, the positive magnetoresistance 
observed below T/v = 88 K in Figure 1(b) can be under- 
stood in terms of the orbital motion of charge carriers 
within this stripe-like ground state. The increased resis- 
tivity for H \\ c may be caused by enhanced scattering 
for charge motion along one-dimensional chains due to 
the Lorentz force. Understanding the interplay between 
spin order and charge motion in this correlated electron 
system remains an important issue for further study. 
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